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1.  Introduction 

This  report  covers  our  activity  during  the  period  of  June  1,1988  through  November  30,1988. 
Progress  has  been  made  along  two  major  lines  of  research,  namely  : 

a)  structure/control  synthesis  of  fiber  composite  lifting  surfaces  in  the  subsonic/ supersonic  flow 
regimes; 

and 

b)  development  of  a  transonic  acroservoelastic  capability. 

Considering  the  multidisciplinary  nature  of  the  research  and  its  wide  scope,  research  activity  in 
the  past  six  months  was  still  directed  at  the  programming,  checking  and  studying  the  selected 
methods  of  analysis  for  each  of  the  disciplines  considered.  At  the  same  time  major  portions  of  the 
programming  needed  to  integrate  the  analysis  techniques  and  behavior  sensitivity  calculations  into 
a  multidisciplinary  optimization  capability  have  been  completed. 

2.  Structure/  Control  Synthesis  of  Lifting  Surfaces  in  the  Subsonic  and  Supersonic  Flow  regimes. 

The  Structural  Model 


Further  testing  of  the  structural  analysis/sensitivity  equivalent  plate  module  has  been  completed. 
After  establishing  its  accuracy  and  computational  efficiency  for  the  analysis  of  practical  wing/ 
control  surface  assemblies,  additional  numerical  testing  was  added  to  verify  its  analytic  behavior 
sensitivity  calculations.  Validity  of  the  analytic  behavior  sensitivities  was  studied  by  comparing 
them  to  finite  difference  derivatives. 

Some  typical  results  are  presented  in  Tables  I  and  2.  Table  1  presents  displacement  and  skin 
quadratic  strength  failure  criterion  derivatives  with  respect  to  skin  thickness  for  the  low  aspect 
ratio  wing  of  Ref.  1.  Displacement  and  spar  cap  stress  sensitivities  for  the  wing  of  Ref.  2  with 
respect  to  skin  thickness  and  to  spar  cap  area  are  shown  in  Table  2.  Basted  on  such  behavior 
sensitivities,  the  robustness  of  alternative  approximations  used  for  different  constraints  can  be 
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assessed.  Fig.  1  is  a  typical  example  of  such  a  study,  where  the  constraint  is  an  equivalent  stress  ( 
calculated  by  a  quadratic  combined  stress  criterion  for  the  skin  of  a  low  aspect  ratio  wing  )  and 
the  design  variable  is  the  skin  thickness.  Among  the  alternative  approximations  (  in  terms  of 
different  intermediate  design  variables  )  the  well  known  superiority  of  the  reciprocal  intermediate 
design  variable  for  this  type  of  constraint  is  evident  from  Fig.  1 . 

Overall,  the  equivalent  plate  computer  structural  modeling  capability  has  proven  to  be  accurate 
and  fast  for  both  analysis  and  behavior  sensitivity  calculations.  A  detailed  documentation  of  the 
structural  capability  has  been  completed  and  a  preliminary  copy  is  submitted  with  this  report 
(Ref.  3). 

Finite  State  Modeling  of  Unsteady  Aerodynamic  Loads. 

The  generalized  aerodynamic  loads  in  the  Laplace  transformed  small  perturbation  equations  of 
motion  are  transcendental  functions  of  the  Laplace  variable  s.  To  use  modem  control  system 
analysis  and  design  techniques,  it  is  necessary  to  cast  them  in  Linear  Time  Invariant  (LTI)  state 
space  form.  The  common  practice  is  to  match  rational  function  approximations  to  generalized 
aerodynamic  loads  calculated  for  harmonic  motion  at  a  set  of  reduced  frequencies  (Ref.  4).  There 
is  a  resulting  increase  in  the  order  of  the  LTI  state  space  model  due  to  the  addition  of 
aerodynamic  states.  This  increase  in  size  can  be  quite  significant  and  make  it  computationally 
expensive  to  carry  out  control  system  analysis  and  behavior  sensitivity  analysis  using  state  space 
techniques.  In  the  Minimum  State  Method  (  MSM  )  (Ref.  5  )  the  form  of  the  approximant  is 
chosen  so  that  the  number  of  added  Mates  needed  to  achieve  a  given  quality  of  fit  is  significantly 
reduced  compared  to  other  method-*  Order  reduction  of  the  state  space  model  used  for 
aeroservoelastic  stability  and  control  analysts  is  essential  for  synthesis  purposes  because  it 
produces  computationally  fast  unuh'is  o*ks  Ihus  the  Minimum  State  Approximation  is  ideally 
suited  for  a  situation  where  optimization  requirements  lead  to  repetitive  calculations,  and  it  was 
selected  for  the  finite  state  unsteady  uvrody  numic  model  employed  in  our  research. 
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However,  it  should  be  emphasized  that  very  limited  experience  exists  on  the  application  of  this 
method  to  practical  airplane  configurations.  In  absence  of  modal  reduction,  which  is  one  of  the 
features  of  our  research,  when  the  flutter  analysis  is  carried  out  with  the  full  order  generalized 
mass,  stiffness  and  aerodynamic  matrices  corresponding  to  the  same  set  of  polynomial  Ritz 
functions  used  for  structural  analysis,  there  is  no  previous  experience  with  this  method.  Therefore 
a  computer  code,  capable  of  implementing  this  method  and  testing  its  accuracy,  had  to  be 
developed.  Results  from  this  code  had  to  be  compared  with  other  available  methods  used  for  the 
finite  state  modeling  of  unsteady  aerodynamic  loads.  The  results  obtained  during  these  tests  were 
excellent. 

For  example,  a  44  x  44  matrix  of  generalized  aerodynamic  forces  for  the  YF16  airplane  with  tip 
missiles  is  approximated  successfully  using  only  22  lag  terms.  Comparison  with  a  one  lag  term 
Roger  approximation  (which  will  add  44  aerodynamic  states  to  the  model  ,  see  Ref.  4)  is  shown 
in  Fig  2.  As  the  results  indicate,  the  one  term  Roger  approximant  is  not  capable  of  capturing  the 
frequency  dependent  complex  behavior  of  generalized  aerodynamic  force  terms,  while  the 
adequacy  of  the  MSM  approximation  is  clearly  demonstrated. 

Classical  Flutter  Analysis  based  on  the  MSVl  Approach. 

To  study  additional  aspects  of  integration  of  the  equivalent  plate  structural  approach,  the  PCKF 
(Piecewise  Continuous  Kernel  Function)  aerodynamics  (see  previous  progress  reports  and  Ref.  6) 
and  the  MSM  finite  state  modeling  method,  classical  flutter  calculations  have  been  carried  out  for 
several  wings.  The  swept  wing  of  Ref  '  is  a  typical  example.  A  ten  term  complete  polynomial 
series  was  used  as  a  Ritz  scries  for  ihc  v.intdcvered  wing.  Generalized  aerodynamic  forces  were 
generated  for  this  set  of  generalized  loordinjtcs  and  fitted  using  the  MSM  approximant  with  four 
added  states.  A  typical  root  locus  u  Town  in  Fig.  3  for  the  case  where  the  model  order  is 
reduced  using  3  natural  modes  of  vibration  It  shows  that  there  is  coupling  between  structural 
and  aerodynamic  poles  with  the  MSM  model,  although  those  aerodynamic  added"  poles  which 
are  complex  are  always  more  heavily  damped  than  the  structural  poles  in  the  same  frequency 
range.  The  behavior  of  poles  near  the  imaginary'  axis  is  predicted  very  well. 
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When  no  modal  reduction  is  carried  out  and  the  root  loci  are  calculated  with  the  full  order 
structural  model,  the  following  lessons  have  been  learned  : 

Since  the  MSM  approximation  is  valid  only  up  to  the  highest  reduced  frequency  used  in  the 
aerodynamic  calculations  its  use  beyond  this  frequency  range  produces  questionable  results.  When 
the  full  order  structural  model  is  used,  without  modal  reduction,  some  very  high  frequencies  are 
included  among  the  structural  poles.  Modified  by  the  aerodynamics  at  these  frequencies,  they  can 
lead  to  unreliable  prediction  of  the  high  frequency  poles.  This  problem  can  be  resolved  in  a 
straight  forward  manner.  High  frequency  poles  are  questionable  because  of  structural  modeling 
limitations.  Thus  satisfactory  aeroscrvoelastic  stability  analysis  is  achieved  by  ignoring  all  complex 
eigenvalues  whose  damped  frequency  is  above  the  range  of  frequencies  used  with  the  MSM 
approximation.  In  a  full  order  flutter  analysis  carried  out  for  the  swept  wing  of  Ref.  7  the 
behavior  of  the  low  frequency  poles  is  essentially  the  same  as  that  predicted  by  the  reduced  order 
model. 

Currently  we  are  studying  the  calculated  pole  structure  for  large  full  order  structural  models  (the 
YF 16  model  with  forty  degrees  of  freedom).  Based  on  the  results  of  this  study  a  decision  will  be 
made  whether  to  use  modal  reduction  or  full  order  models  for  the  aeroservoelastic  analysis.  Both 
options  are  included  in  the  computer  program  which  has  been  developed. 

3.  Development  of  Transonic  Acroelastic  Capability. 

A  transonic  aeroservoelastic  capability  is  currently  in  an  advanced  stage  of  programming  and 
development.  The  program  is  based  on  a  combination  of  the  LISSA  structural  analysis  code 
developed  in  the  course  of  this  research  and  a  modified  version  of  the  XTRAN2L  code  which  is 
being  extensively  modified  so  as  to  include  unsteady  three  dimensional  finite  span  corrections. 

It  is  important  to  emphasize  that  the  original  XTRAN2L  code  described  in  Ref.  8  is  not  an 
aerodynamic  code;  it  is  more  of  an  acroelastic  code,  since  the  structural  equations  and  the 
aerodynamic  equations  arc  integrated  simultaneously.  For  our  purposes  only  the  aerodynamic 
pressure  calculations  in  the  time  domain  are  required.  Thus  all  nonessential  portions  of  the  code 
such  as  forced  harmonic  response  and  pitch  pulse  routines  have  been  deleted,  and  a  three 
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dimensional  correction  to  the  program  has  been  implemented.  The  3-D  corrections  are  based  on  a 
modified  version  of  an  inverse  design  approach.  In  order  to  have  a  continuous  unsteady  pressure 
coefficient  distribution  over  the  wing  surface,  a  surface  spline  technique,  described  in  Ref.  9,  with 
a  scaled  spanwise  coordinate  is  used.  The  program  also  utilizes  on  line  parameter  estimation  to 
determine  flutter  boundaries. 

Numerous  portions  of  the  analysis,  consisting  of  program  subroutines,  have  already  been 
developed  and  tested.  Presently  these  elements  are  being  combined  with  each  other  and  with  the 
LISSA  structural  analysis  code.  In  this  combination  the  structural  and  aerodynamic  grid  points 
have  to  be  treated  carefully  to  ensure  appropriate  compatibility. 

4.  Status 


Table  3  shows  the  current  status  of  research  activities  associated  with  the  development  of  the 
analysis  and  sensitivity  capabilities  for  the  multidisciplinary  synthesis  of  wings  for 
subsonic/supersonic  flight.  A  paper  titled  :  'An  Integrated  Approach  to  the  Optimum  Design  of 
Actively  Controlled  Composite  Wings"  was  presented  at  the  Second  NASA/  Air  Force 
Symposium  on  Recent  Advances  in  Multidisciplinary  Analysis  and  Optimization,  held  at 
Hampton,  Virginia,  September  28-80,1988.  The  paper,  which  is  in  a  Viewgraph  /  Notes  format, 
constitutes  Appendix  A  of  this  progress  report.  A  second  more  comprehensive  paper  is  being 
prepared  for  presentation  at  the  forthcoming  AIAA/ASME/ASCE/AHS/ASC  30th  Structures, 
Structural  Dynamics  and  Materials  Conference  in  Mobile,  Alabama.  We  expect  the  transonic 
aeroservoelastic  code  to  become  operational  around  the  middle  of  January  1989. 
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FIGURE  1  : 


APPROXIMATION  OF  A  QUADRATIC  STRESS  CRITERION  AT  A  POINT  ON  THE 
SKIN  OF  A  LOW  ASPECT  RATIO  WING  (REF.  1) 


FIGURE  2 


SOME  RECENT  EXAMPLES  OF  QUALITY  OF  FIT  FOR  A  YF16  COMPLETE  A/C 
CONFIGURATION  : 

(  44  POLYNOMIAL  GENERALIZED  COORDINATES  ) 

SOME  LOW  ORDER  SHAPE  FUNCTIONS  FOR  THE  WING  BOX  : 

A  (X^)«  I.  A  (jcjr)-  x,  A  (*»?)-  X* 

A  HIGHER  ORDER  SHAPE  FUNCTION  :ftt  (xjO-  x* 

ROGER  APPROXIMATION  BASED  ON  1  LAG  (44  ADDED  STATES) 

MINIMUM  STATE  BASED  ON  22  LAGS  (22  ADDED  STATES) 


o.oso  ■ 


IMAG  A(l,4) 


IMAG  A<U1) 


IMAG  A(U) 


ROGER  - 

MIN  STATE 
•EXACT  \ 


0.025 


0.000 


0.0012 


0.0000 


l  ROGER  | 
j  MIN  STAT 
N EXACT 


-0.023 


-0.0012 


REAL  A(  1 ,2) 

.  1  1"~  . . 

_n  «  -0.6  -0.4  -0.2 


-0.050  ■ 


REAL  A(  1.4) 

H - - - -T- - - - 

-0.25  -0.20  -0.15  -0.10  - 


-0.0023 


-0.010 


REAL  A(  1.11) 


-0.003 


■  iGL'RE  3  : 


i  *)  3  natural  modes  in  the  range  of  frequencies  corresponding  to  valid  Min  State  approximation. 


table  I 


Analytical  versus  Finite  Difference  Sensitivities  of  the  Giles’  Wing  Deflections  and  Stresses. 


Behavior 

Vnalytic 

Finite 

Finite 

Function 

Dcnvativet  ♦) 

Difference 

Derivative!  •) 

(A  =  1.0  %1 

Difference 

Derivative!  * ) 

(A  -  0.1  %  1 

Tip  Displacement 

Constraint 

at  x*  1 50'. y  =  200' 

■24  0474 

*23.8362 

-24.C230 

Quadratic  Stress 

Constraint 

at  x  =  M"  y  =  40' 

•4  8814 

-9.7876 

-9  8710 

(*)  Design  Variable  :  t  .the  thickness  of  cover  skin  of  the  outer  wing  section. 
5  psi  load.  4  '  max  displacement,  aluminum  skins. 
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TABLE  2  : 

Analytical  versus  finite  Difference  Sensitivities  01  the  Denkc  Wing  i  Displacements  and  Stresses i. 


Design  Variable  .  Skin  Thickness  It) 


Behavior 

Function 

Analytic 

Derivative 

Finite 

Difference 

Derivative 

(A  -  1.0  %) 

Finite 
Difference 
Derivative 
(A -0.1  %  ) 

Rear  Spar 

Tip 

Displacement 

-14.S476 

-14.5294 

-14.5500 

Rear  Spar  Cap  -15.8700 

Stress  at  the 

Root  ( normalized  i 

- 15.5339 

-15.6329 

Design  Variable  : 

Rear  Spar  Cao  Area  tin1) 

Behavior 

Function 

Analytic 

Derivative 

Finite 

Difference 

Derivative 

(A  -  1.0  %) 

Finite 

Difference 

Derivative 

(A -0.1  %) 

Rear  Spar 

Tip 

Displacement 

- 1  2314 

-1.2272 

-1.2318 

Rear  Soar  Cap 

- 1  2796 

-1.2753 

-1.2791 

Stress  at  the 
Root  (normalized i 


TABLE  3  : 


ANALYSIS  AND  BEHAVIOR  SENSITIVITY 
STATUS 


STRUCTURE 


FORMULATION  :  ANALYSIS  + 

SENSITIVITY  + 

ANALYSIS  IMPLEMENTATION  :  + 

ANALYSIS  TESTING  :  + 

SENSITIVITY  IMPLEMENTATION  :  + 

SENSITIVITY  TESTING  :  + 

APPROXIMATION  CONCEPTS 
PERFORMANCE  ASSESSED  :  + 


AERODYNAMICS  * 

CONTROL 

+■ 

+ 

+ 

4- 

+■ 

+ 

in  progress 

in  progress 

in  progress 

AERODYNAMICS  INCLUDE  : 

*  UNSTEADY  AERODYNAMICS  FOR  SERVOAEROELASTIC  ANALYSIS 

•  STEADY  TRIM  AND  DRAG  CALCULATIONS 


APPENDIX  A 
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INTRODUCTION  AND  OUTLINE  OF  THE  PRESENTATION 


The  importance  of  interactions  among  the  various  disciplines  in  airplane  wing  design  has  been 
recognized  for  quite  some  time.  With  the  introduction  of  high  gain,  high  authority  control  systems 
and  the  design  of  thin,  flexible,  light  weight  composite  wings,  the  integrated  treatment  of  control 
systems,  flight  mechanics  and  dynamic  acroelasticity  became  a  necessity.  A  research  program  is 
underway  now  aimed  at  extending  structural  synthesis  (Ref.  1)  concepts  and  methods  to  the 
integrated  synthesis  of  lifting  surfaces,  spanning  the  disciplines  of  structures,  aerodynamics  and 
control  for  both  analysis  and  design.  Mathematical  modeling  techniques  are  carefully  selected  to 
be  accurate  enough  for  preliminary  design  purposes  of  the  'complicated,  built  up  lifting  surfaces 
of  real  aircraft  with  their  multiple  design  criteria  and  tight  constraints'  (Ref.  2,p.l7).  The 
presentation  opens  with  some  observations  on  the  multidisciplinary  nature  of  wing  design.  A  brief 
review  of  some  available  state  of  the  art  practical  wing  optimization  programs  and  a  brief  review 
of  current  research  effort  in  the  field  sen  e  to  illuminate  the  motivation  and  support  the  direction 
taken  in  our  re  search. (These  reviews  are  not  exhaustive,  and  the  interested  reader  is  referred  to  the 
review  papers  Refs.  3-8.)  The  goals  of  this  research  effort  will  be  presented  next,  followed  by  a 
description  of  the  analysis  and  behavior  sensitivity  techniques  used.  The  presentation  will 
conclude  with  a  status  report  and  some  forecast  of  upcoming  progress. 


OUTLINE  OF  THE  PRESENTATION 


BRIEF  REVIEW  OF  CURRENT  WING  OPTIMIZATION  CAPABILITIES 
AND  RESEARCH  ACTIVITY,  SOME  OBSERVATIONS 


GOALS  FOR  MULTIDISCIPLINARY  WING  SYNTHESIS  RESEARCH  AT  UCLA 


DESCRIPTION  OF  ANALYSIS  TECHNIQUES  CHOSEN 


STATUS  REPORT  ON  THE  SYNTHESIS  CAPABILITY  UNDER  DEVELOPMENT 


THE  MULTIDISCIPLINARY  NATURE  OF  WING  DESIGN 


Figure  1  describes  the  multidisciplinary  nature  of  wing  design.  Discussion  is  limited  to  wings 
operating  in  the  subsonic  to  low  supersonic  flight  speeds,  so  that  thermal  effects  can  be  neglected. 
It  is  instructive  to  unite  the  sets  of  Preassigned  Parameters  and  Design  Variables  (Ref.l)  into  the 
set  of  "Design  Parameters',  whose  elements  define  a  particular  wing  design.  Which  of  the 
parameters  will  be  preassigned  and  which  will  be  used  as  design  variables  depends  on  the  level  of 
application  for  optimization  techniques  in  the  hierarchy  described  in  Ref.  1,  namely,  whether  the 
design  space  includes  sizing,  configuration(geometry)  or  topological  design  variables.  The  set  of 
behavior  functions,  from  which  constraints  and  objectives  will  be  selected,  can  be  divided  into  two 
categories.  Primary  (system  level)  Behavior  Functions  are  those  performance  measures  which 
determine  the  overall  quality  and  competitiveness  of  the  wing.  Secondary  (sub-system  level) 
Behavior  Functions  are  the  behavior  functions  which  must  be  taken  into  account  during  the 
design  to  guarantee  the  prevention  of  failure  in  all  possible  failure  modes  and  introduce  known 
constraints  on  subsystem  performance.  They  are  usually  not  the  real  design  objectives  although 
sometimes  there  is  high  correlation  between  a  secondary  behavior  and  a  primary  behavior 
function  (e.g.  mass  and  airplane  performance). 


\ 
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THE  MULTIDISCIPLINARY  NATURE  OF  WING  DESIGN 


DISCIPLINE 


CONFIGURATION 

DESIGN 

PARAMETERS 


TOPOLOGICAL 
DESIGN 
PARAMETERS  : 


SECONDARY 
BEHAVIOR 
FUNCTIONS  : 


PRIMARY 

BEHAVIOR 

FUNCTIONS: 


STRUCTURES/ 

STRUCTURAL 

DYNAMICS 


PLAN  FORM  SHAPE 
(SWEEP.AR. 

TAPER  RATIO) 
AIRFOIL  CROSS 
SECTION,  PLY 
ANGLE 


NUMBER  OF  RIBS, 
SPARS,  PLIES 


AERODYNAMICS 


CAMBER, 

LE/TE  CONTROL 

SURFACE 

DEFLECTION 


PLANFORM  SHAPE 
(SWEEP, AR 
TAPER  RATIO) 
AIRFOIL  CROSS 
SECTION 


NUMBER  OF 
LE,  TE  CONTROL 
DEVICES 


DEFLECTION/SLOPE, 
STRESS. BUCK  LING, 
NATURAL  FREQ'S, 
FATIGUE. 

LOAD  FACTOR. 
MASS, 

CONTROL 

EFFECTIVENESS 


DRAG, 

LIFT.CLMAX, 
DRAG  POLAR, 


CONTROL 

EFFECTIVENESS 


.  AIRPLANE 
PERFORMANCE  : 
(ROLL  RATE,  TURN  RATE, 
ENDURANCE,  RANGE, 
COST...  ) 


SIZING 

DESIGN  PARAMETERS  : 

(DESIGN 

VARIABLES  OR  ELEMENT  SIZE 
PREASSIGNED)  (AREA, 
PARAMETERS  THICKNESS) 


Figure  2 


CONTROL 


GAINS, 

TRANSFER 

FUNCTION 

COEFFICIENTS 


ORDER  OF 
TRANSFER 
FUNCTIONS 


CONTROL 

SYSTEM 

STRUCTURE 

AND 

CONNECTIVITY 


DYNAMIC 

STABILITY 

(INCLUDING 

FLUTTER, 

BODY  FREEDOM 
FLUTTER), 

HINGE  MOMENTS, 

CONTROL 

ACTIVITY 


-3- 


SOME  EXISTING  PROGRAMS  FOR  PRACTICAL  WING  OPTIMIZATION 


Several  approaches,  with  a  varying  degree  of  multidisciplinary  capability,  aimed  at  the  synthesis  of 
practical  composite  wings  were  developed  during  the  seventies  (Ref.  3-8).  In  addition  to  the 
constraints  on  stress,  displacement  and  acroelastic  stability,  performance  constraints  in  terms  of 
induced  drag  or  drag  polar  specification  were  added  in  the  TSO  computer  code  (Ref.  9,10)  and  to 
WIDOW  AC  (Ref.  12-14).  It  was  reported  recently  that  a  rudimentary  scrvoaeroelastic  analysis 
capability  was  about  to  be  inserted  into  the  ASTROS  computer  code  (ref.  20).  It  should  be 
noticed  that  except  for  the  TSO  code,  the  design  space  in  the  programs  contain  only  structural 
design  variables,  thus  they  are  really  multidisciplinary  in  analysis  only.  The  TSO  code  makes  it 
possible  to  include  some  configuration  design  variables  (the  fiber  orientation  of  cover  skin  layers) 
and  some  aerodynamic  constraints  in  the  form  of  wing  twist  or  camber  distribution  under  load. 
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SOME  EXISTING  PROGRAMS  FOR  PRACTICAL  WING  OPTIMIZATION 

PROGRAMS  SURVEYED  : 

♦  TSO.  FASTOP,  WIDOWAC,  ELFIM,  ASTROS 

ANALYSIS  PROBLEM 
DISCIPLINES  : 

♦  STRUCTURES,  AERODYNAMICS,  AEROELASTICITY. 

MODELING  : 

STRUCTURAL : 

♦  EQUIVALENT  PLATE  (TSO) 

♦  FINITE  ELEMENTS  (FASTOP, WIDOWAC, ELFINI,ASTROS) 

AERODYNAMIC 
STEADY  AERODYNAMICS  : 

*  LINEAR  POTENTIAL  PANEL  METHODS  (e  g.  Woodward  in  TSO) 
UNSTEADY  AERODYNAMICS  : 

*  DOUBLET  LATTICE  (  M  <  1  )  (TSO,  FASTOP,  ASTROS) 

♦  KERNEL  FUNCTION  (  M  <  1  )  (WIDOWAC) 

*  MACHBOX  (  M  >  I  )  (FASTOP) 

♦  POTENTIAL  GRADIENT  METHOD  (  M  >  1  )  (ASTROS) 

♦  PISTON  THEORY  (  M>  >  I  )  (WIDOWAC) 

BEHAVIOR  SENSITIVITY 

♦  ANALYTIC  (FASTOP,  WIDOWAC,  ELFINI,  ASTROS) 

♦  FINITE  DIFFERENCES  (TSO) 

SYNTHESIS  PROBLEM 
PREASSIGNED  PARAMETERS  : 

♦  PLANFORM  SHAPE,  CROSS  SECTION,  STRUCTURAL  TOPOLOGY,  MATERIALS 
DESIGN  VARIABLES  : 

♦  STRUCTURAL  SIZING  (TSO  -  ALSO  FIBER  ORIENTATION  ) 

CONSTRAINTS  : 

♦  DEFLECTION,  STRESS,  FLUTTER,  DIVERGENCE,  CONTROL  EFFECTIVENESS 

♦  DRAG  (TSO,  WIDOWAC) 

♦  BUCKLING  (ELFINI,  ADDED  TO  WIDOWAC  IN  A  SIMPLIFIED  FORM) 

OBJECTIVE  FUNCTIONS  : 

♦  MASS,  DRAG,  CONTROL  EFFECTIVENESS 

OPTIMIZATION 

♦  MATH  PROGRAMMING  (TSO) 

♦  MATH  PROGRAMMING  +  APPROXIMATION  CONCEPTS  (WIDOWAC,  ELFINI, 
ASTROS) 

♦  OPTIMALITY  CRITERIA  (FASTOP) 


THE  NEED  FOR  MULTIDISCIPLINARY  WING  OPTIMIZATION 


During  the  last  decade  structural  synthesis  has  matured.  Realistic  designs  described  by  a  large 
number  of  design  variables  and  subject  to  a  variety  of  load  conditions  can  now  be  efficiently 
treated.  However,  it  is  still  quite  common  to  find  fixes  and  modifications  being  introduced  late  in 
the  development  stage  of  fighter  aircraft,  when  aeroservoelastic  effects,  rigid  body-  elastic  mode 
coupling  or  static  aeroelastic  effects  have  not  been  properly  accounted  for  in  the  design  process 
(Ref.  21-27).  At  the  same  time,  following  almost  twenty  years  of  progress  in  active  flutter 
suppression  and  gust  alleviation  (Ref.  28-32),  there  is  a  growing  recognition  that  multidisciplinary 
interactions  might  be  harnessed  to  benefit  modem,  complex  wing  designs.  However,  a  review  of 
the  literature  reveals  that  the  application  of  modem  optimization  methods  to  wing  design 
problems  involving  multiple  objective  functions  and  a  diverse  mix  of  constraints  based  on 
analyses  from  several  discipline  areas  (e.g.  structures,  structural  dynamics,  controls,  aerodynamics 
and  performance)  has  not  yet  been  treated  in  a  comprehensive  and  realistic  manner.  To 
overcome  the  inherent  complexity  and  address  the  computationally  intensive  nature  of  this 
problem  two  approaches  have  been  suggested  in  the  literature.  The  first  approach  is  based  on  the 
application  of  multi-level  decomposition  techniques  combined  with  existing  tools  for  detailed 
analysis  and  sensitivity  analysis  for  each  of  the  disciplines  (Ref.  40,41).  The  second  approach 
attempts  to  gain  some  insight  into  the  nature  of  the  problem  by  using  highly  simplified 
mathematical  models  or  simple  airplane  configurations  for  structural,  aerodynamic  and  control 
system  analysis  (Ref.33-39).  Research  is  now  underway  in  several  research  centers  and  universities 
in  two  main  directions  : 

a)  the  addition  of  control  system  sizing  type  design  variables  to  a  design  space  spanning  design 
variables  for  structures  and  control  (control  augmented  aeroelastic  optimization)  (Ref.  35,42) 

b)  expanding  the  wing  design  space  by  adding  configuration  design  variables  (structural  and 
aerodynamic  shape)  (Ref.  37-39). 
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THE  NEED  FOR  MULTIDISCIPLINARY  WING  OPTIMIZATION 


MOTIVATION  : 

PROBLEMS  : 

*  SERVOAERELASTIC  INTERACTIONS 

*  INSUFFICIENT  CONTROL  EFFECTIVENESS 

*  BODY  FREEDOM  FLUTTER 

HARNESS  MULTIDISCIPLINARY  INTERACTIONS  FOR  ACHIEVING  BETTER 
DESIGNS 


CURRENT  DIRECTIONS  IN  RESEARCH  : 

OPENING  UP  A  SIZING  DESIGN  SPACE  TO  INCLUDE  SIZING  TYPE  CONTROL 
SYSTEM  DESIGN  VARIABLES 

OPENING  UP  THE  DESIGN  SPACE  TO  INCLUDE  STRUCTURAL  AND 
AERODYNAMIC  SHAPE  DESIGN  VARIABLES 


PROBLEMS  : 

♦  HEAVY  COMPUTATIONAL  COST  OF  ANALYSIS 

♦  LACK  OF  INTUITION  AND  EXPERIENCE  TO  GUIDE  IN  CONSTRUCTING 
ROBUST  APPROXIMATIONS 

♦  ABSENCE  OF  STANDARD  TERMINOLOGY,  CRITERIA,  MODELING  AND 
ANALYSIS  METHODS  WIDELY  ACCEPTED  IN  ALL  DISCIPLINES 


APPROACHES  TO  THE  PROBLEM  : 

*  SIMPLIFIED  MODELING  (BEAM, STRIP  THEORY)  OR  STUDIES  INVOLVING 
SIMPLE  CONFIGURATIONS  (SAILPLANES)  - 

ONE  LEVEL  OPTIMIZATION 
GAIN  INSIGHT 

*  MULTILEVEL  DECOMPOSITION  BASED  ON  DETAILED  MODELING  AND 
SENSITIVITY  ANALYSIS 


Figure  4 
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RESEARCH  GOALS 


In  Ref.  2  Ashley  w  rites  :*  In  the  absence  of  experience  w  hen  new  technology  is  being  tried  for  the 
first  time,  the  search  for  extremals  can  produce  unanticipated,  surprising  and  often  very 
satisfactory  discoveries".  But  he  adds  a  word  of  caution  :  'Yet  the  counterintuitive  may  also  be 
counterproductive  and  even  ridiculous.  Very  undesirable  consequences  can  result  from  omission 
or  careless  handling  of  constraints". 

It  is  one  of  the  major  goals  of  the  present  research  to  begin  to  bridge  the  gap  between  over 
idealized  modeling  and  detailed  structural  and  aerodynamic  modeling  by  introducing  balanced 
design  and  analysis  models  that  capture  the  essential  behavior  characteristics,  without  making  the 
integrated  multidisciplinary  design  optimization  task  intractable.  This  balanced  approach 
combines  high  quality,  approximate,  but  computationally  efficient  analyses  for  the  structural, 
aerodynamic  and  aeroservoelastic  behavior  of  realistic  composite  wings.  Thus,  the  entire 
optimization  problem  may  be  treated  at  one  level  without  the  need  for  multilevel  decomposition. 
A  rich  variety  of  constraints  makes  it  possible  to  study  the  efTect  of  multidisciplinary  interactions 
on  synthesis  as  well  as  on  analysis. 
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R ESEARCH  OBJECTIVES  AND  SELECTED  APPROACH 


OBJECTIVES  . 


DEVELOP  MULTIDISCIPLINARY  WING  SYNTHESIS  CAPABILITY  WITH  AN 
EMPHASIS  ON  STRUCTL'RE/CONTROL/UNSTEADY  AERODYNAMICS 
INTERACTION 

BRIDGE  THE  GAP  IN  MODELING  DETAIL  BETWEEN  THE  VERY  SIMPLE  AND 
DETAILED  ANALYSIS  TECHNIQUES  SO  AS  TO  ENABLE  MULTIDISCIPLINARY 
SYNTHESIS  OF  REAL  WINGS  FOR  PRELIMINARY  DESIGN 

STUDY  THE  CONSTRUCTION  OF  ROBUST  APPROXIMATIONS  TO  BEHAVIOR 
FUNCTIONS 

PROVIDE  A  TEST  CASE  FOR  ASSESING  DECOMPOSITION  TECHNIQUES 


SELECTED  APPROACH 


ANALYSIS  : 

CAREFUL  SELECTION  OF  ANALYSIS  TECHNIQUES  - 

GOOD  ACCURACY 

HIGH  COMPUTATIONAL  SPEED 

BALANCED  APPROACH 


BEHAVIOR  SENSITIVITY 
ANALYTIC 

SYNTHESIS  PROBLEM  : 

SIZING  : 

STRUCTURAL.AERODYNAMIC  AND  CONTROL  D.V.'s  PLUS  0,  s 
PREASSIGNED  :  SHAPE,  TOPOLOGY 

CONSTRAINTS  :  S  ,  <r,  SERVOAEROELASIC  STABILITY,  CONTROL  POWER 
ALTERNATIVE  OBJECTIVE  FUNCTIONS  :  MASS,  PERFORMANCE  MEASURES 


OPTIMIZATION  STRATEGY 

MATH  PROGRAMMING  +  APPROXIMATION  CONCEPTS  - 
(  10  ANALYSES  PER  OPTIMIZATION  -  GOAL) 


Figure  5 
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ANALYSIS  METHODS  SELECTED 


The  integrated  optimum  design  capability  outlined  here  is  based  on  approximate  analysis 
techniques  for  the  required  disciplines,  which  are  consistent  with  each  other  in  terms  of  accuracy 
and  efficiency  and  lead  to  a  balanced  treatment.  In  the  structures  area,  an  equivalent  plate 
analysis,  as  incorporated  in  the  TSO  computer  code  (Ref.  10)  and  further  generalized  by  Giles 
(Ref.43-45),  is  used.  Although  the  equivalent  plate  approach  for  structural  modeling  of  low 
aspect  ratio  wings  has  been  known  for  many  years,  it  was  Giles  who  recently  showed  that,  <ising 
present  day  computers,  a  single  high  order  power  series  can  be  used  for  approximating 
displacements  over  wing  planforms  made  of  several  trapezoidal  segments  to  obtain  accurate  stress 
as  well  as  displacement  information.  The  simplicity  of  manipulating  simple  power  series  leads  to 
analytic  rather  than  numerical  integration  for  the  mass  and  stiffness  expressions.  With  the  careful 
organization  of  computer  storage  space  and  ordering  of  calculations,  major  savings  can  be 
achieved  in  terms  of  computation  times  and  core  storage  requirements.  The  extended  equivalent 
plate  approach  is  integrated  with  the  PCKFM  (Piecewise  Continuous  Kernel  Function  Method) 
of  Nissim  and  Lottati  for  lifting  surface  unsteady  aerodynamics  (Ref.  48-51).  This  method 
combines  the  power  of  the  doublet  lattice  method  in  dealing  with  pressure  singularities  with  the 
accuracy  and  speed  of  the  kernel  function  method.  Extensive  numerical  experimentation  has 
demonstrated  (Ref.  48)  that  the  PCKFM  method  is  highly  accurate  and  converges  rapidly.  For 
configurations  involving  control  surfaces,  it  is  faster  and  considerably  more  accurate  than  the 
doublet  lattice  method.  Thus,  it  is  especially  suited  for  calculating  the  generalized  unsteady  air 
loads  (on  lifting  surfaces  made  up  of  wing  and  control  surface  elements)  that  are  needed  for  active 
flutter  suppression  and  gust  alleviation  studies. 

For  the  finite  state  modeling  of  the  unsteady  air  loads,  the  Minimum  State  Method  of  Karpel 
(Ref.  52)  is  used  to  generate  accurate  approximations  to  unsteady  generalized  aerodynamic  forces 
with  addition  of  only  a  small  number  of  augmented  states  to  the  mathematical  model  of  the 
aeroservoelastic  system.  In  comparison  with  other  finite  state  modeling  techniques,  the  number  of 
added  states  needed  in  the  minimum  state  method  can  be  smaller  for  the  same  overall  accuracy  of 
approximation  (Ref.  53).  This  leads  to  a  state  space  model  of  lower  order,  thus  reducing  memory 
requirements  and  computation  times  considerably.  The  integrated  servoaeroelastic  system  is 
modeled  as  a  Linear  Time  Invariant  (LTI)  system  and  its  stability  is  examined  by  computing  the 
eigenvalues  of  a  generalized  eigenvalue  problem. 
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ANALYSIS  METHODS  SELECTED 

STRUCTURE . 

EQUIVALENT  PLATE 

AERODYNAMICS  : 

SUBSONIC/SUPERSONIC  LIFTING  SURFACE  PIECEWISE  CONTINUOUS  KERNEL 
FUNCTION  METHOD  (PCKFM)  (NISSIM/LOTTATI) 

UNSTEADY  AERODYNAMIC  FINITE  STATE  MODELING  : 

MINIMUM  STATE  APPROXIMATION 

CONTROL : 

STATE  SPACE  LTI  SYSTEM  MODELING 

Figure  6 
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EQUIVALENT  PLATE  M0DPLL1NG  OF  AIRPLANE/  WING  CONTROL  SURFACE 
ASSEMBLIES  BY  THE  PRESENT  CAPABILITY 


Figure  7  shows  an  airplane  modeled  as  an  assembly  of  flexible  lifting  surfaces.  Each  lifting  surface 
is  modeled  as  an  equivalent  plate  whose  stiffness  is  controlled  by  contribution  from  thin  cover 
skins  (fiber  composite  laminates)  and  the  internal  structure  (spar  and  rib  caps).  Plate  sections  are 
connected  to  each  other  via  stiff  springs  (to  impose  displacement  compatibility  at  attach  points) 
and  flexible  springs  (representing  the  stiffness  of  actuators  and  their  backup  structure).  Each  wing 
section  can  be  made  of  several  trapezoidal  parts  continuously  connected  to  each  other. 
Concentrated  masses  arc  used  to  model  nonstructural  items  and  balance  masses. 

The  present  equivalent  plate  modeling  capability  makes  it  possible  to  efficiently  analyze  combined 
wing  box/control  surface  configurations.  A  wing  assembly  and  a  canard  or  horizontal  tail  may  be 
attached  to  a  fuselage  (modeled  as  a  flexible  beam  or  a  flexible  plate)  to  simulate  complete 
airplane  configurations.  Modeling  detail  of  all  plate  sections  can  be  identical.  Thus  the  degree  of 
detail  in  modeling  control  surfaces  for  analysis  and  synthesis  is  not  limited,  as  is  the  case  in  the 
TSO  code. 
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SOME  ANALYTICAL  ASPECTS  OF  THE  EQUIVALENT  PLATE  APPROACH 


It  is  a  well  known  fact  in  the  numerical  solution  of  partial  differential  equations,  that  the  use  of  a 
simple  polynomial  series  to  approximate  the  solution  in  a  Ritz  or  Galerkin  analysis  leads  to  ill 
conditioning  of  the  problem  matrices  when  it  is  of  an  order  higher  than  a  certain  degree. 
However,  Giles  (Ref.  43,44)  has  shown  that  when  a  simple  polynomial  series  is  used  in  a  Ritz 
solution  of  anisotropic  plate  static  and  dynamic  problems,  accurate  displacements,  stresses  and 
natural  frequencies  can  be  obtained  for  practical  wings  before  ill  conditioning  appear.  His  results 
were  obtained  on  a  CDC  Cyber  173  (60  bit  words).  Our  results  obtained  on  an  IBM  3090 
computer  in  extended  precision  and  on  a  SUN  3/280  computer  using  double  precision  support  his 
findings.  When  the  depth  of  the  wing  and  the  thickness  distribution  of  skin  layers  are  also 
expressed  as  power  series,  it  can  be  shown  that  the  stiffness  and  mass  matrices  are  expressed  as 
linear  combinations  of  certain  area  and  line  integrals  and  polynomial  terms  calculated  at  points 
where  wing  section  are  connected  or  where  concentrated  masses  are  placed.  These  integrals  and 
polynomial  tables  are  fixed  once  a  planform  shape  is  given.  Thus  they  are  evaluated  only  once  at 
the  beginning  of  an  optimization  task.  This  leads  to  major  computation  time  savings  along  with 
the  fact  that  the  relatively  small  number  of  generalized  coordinates  needed  to  accurately 
approximate  displacement  and  stresses  in  a  wing  section  (about  21-30)  result  in  small  mass  and 
stiffness  matrices! although  fully  populated)  compared  with  finite  element  analysis. 
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SOME  ANALYTICAL  ASPECTS  OF  THE  EQUIVALENT  PLATE  APPROACH 


POLYNOMIAL  FUNCTIONS  : 

THE  SHAPE  FUNCTIONS  :  fix j/)  =  x*  y  (*  =  "<«).  n  =  *0 ) 


A  TYPICAL  SKIN  LAYER  THICKNESS  DISTRIBUTION  :  t(xy)=ZT.  x *  / 
(kj  DEPEND  ON  i) 

WING  DEPTH  :  Hxy)  =  £//,  x*  y  (r.s  DEPEND  ON  i) 

l»l 

SPAR/RIB  CROSS  SECTIONAL  AREA  :  A(x)-Au+Axx 


FUNDAMENTAL  INTEGRALS  : 

AREA  INTEGRAL  OVER  A  SKIN  TRAPEZOIDAL  SECTION  :  /.- / /  x"  y  dxdy 
LINE  INTEGRAL  OVER  THE  LENGTH  OF  A  RIB  :  /**-  /m  x-  dx 

J  *f 

LINE  INTEGRAL  OVER  THE  LENGTH  OF  A  SPAR  :  ls~  =  /£»*  x(y>*  y  dy 
ASSEMBLY  : 

MASS  AND  STIFFNESS  MATRICES  ARE  LINEAR  COMBINATIONS  OF  THESE 
INTEGRAL  TERMS  WITH  VARYING  INDICES  :  (m,n) 


DISPLACEMENT  RITZ  SERIES  FOR  TWO  WING  SECTIONS  CONNECTED  TO  EACH 
OTHER  VIA  SPRINGS  : 

wi  =  [ . (* ly-'  W . ]  {?>}  =  [ . (x2)"»  {yTf*  .]  {*} 

FUNDAMENTAL  POLYNOMIAL  TERM  TABLES  : 


m  n  m  n  m  n  m  n 

*\y\  .  x2y2  ,  xty2  ,  x2yx 
ASSEMBLY  : 

ATTACHMENT  CONTRIBUTION  TO  STIFFNESS  MATRIX  IS  A  LINEAR 
COMBINATION  OF  POLYNOMIAL  TERMS  TAKEN  FROM  THE  FUNDAMENTAL 
TABLES. 


Figure  8 
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FEATURES  OF  THE  PRESENT  EQUIVALENT  PLATE  MODELLING 


In  order  to  structurally  analyze  (statics  and  dynamics)  wing/control  surface/  canard  or  tail 
configurations  and  to  accelerate  the  generation  of  approximate  problems  for  synthesis,  the 
equivalent  plate  approach  of  Giles  was  further  extended  to  include  multi-element  wing 
box/control  surfaces  plus  analytic  behavior  sensitivity  derivatives  with  respect  to  structural  design 
variables.  Stiffness  and  mass  matrices  can  now  be  generated  using  analytic  integration  for  wing 
structures  made  of  composite  skins,  spars  and  ribs,  concentrated  masses  and  equivalent  springs 
which  connect  plate  sections  to  each  other. 
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FEATURES  OF  THE  PRESENT  EQUIVALENT  PLATE  MODELLING 


CONFIGURATIONS  MODELLED  INCLUDE  : 

WING/  CONTROL  SURFACE/  CANARD/  FUSELAGE  ASSEMBLIES 
FUSELAGE  AND  MISSILES  CAN  BE  MODELLED  AS  EQUIVALENT  BEAMS 
DESIGN  VARIABLES  INCLUDE  : 

SKIN  LAYER  THICKNESS  DISTRIBUTION  POLYNOMIAL  COEFFICIENTS, 

SPAR /  RIB  CAP  AREA  DISTRIBUTION  (LINEAR  ALONG  SPAR/RIB  LINE) 
CONCENTRATED  MASSES 

LINEAR  AND  ROTATIONAL  SPRING  STIFFNESSES 

ANALYSIS  CAPABILITY  : 

FAST  STIFFNESS, MASS  MATRIX  GENERATION 

STATIC  SOLUTION  FOR  DISPLACEMENT!:  AND  STRESSES  UNDER  GIVEN  LOADS 
CALCULATION  OF  NATURAL  FREQUENCIES  AND  MODE  SHAPES 

SENSITIVITY  : 

ANALYTIC  BEHAVIOR  SENSITIVITY  ANALYSIS  FOR  DISPLACEMENTS,  SLOPE, 
QUADRATIC  FAILURE  CRITERIA  FOR  STRESSES  IN  SKINS.  STRESSES  IN 
SPAR/RIB  CAPS 

ADJOINT  OR  DIRECT  METHOD  -  OPTIONAL 


Figure  9 
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NUMERICAL  TESTING 


Extensive  numerical  tests  were  carried  out  to  study  the  accuracy  of  the  present  equivalent  plate 
modeling  and  asses  its  computational  efficiency.  Several  wings  of  different  construction,  aspect 
ratio  and  thickness  were  used.  Displacements,  stresses  in  skins  and  spar  caps  as  well  as  natural 
frequencies  and  mode  shapes  were  compared  to  finite  element  results  and  to  test  results  where 
available.  As  an  example,  Figure  10  includes  a  comparison  between  YF16  wing  natural 
frequencies  calculated  using  a  detailed  finite  element  analysis,  the  TSO  program  and  our  present 
structural  module.  The  YF16  wing  configuration  includes  a  wing  box  plus  a  leading  edge  flap  and 
a  flaperon.  The  results  demonstrate  the  accuracy  of  the  new  multi-element  equivalent  plate 
modeling  capability  in  analyzing  wing /  control  surface  configurations.  Some  ground  vibration  test 
results  available  in  Ref.  54  made  it  possible  to  check  the  accuracy  of  the  present  code  when  a 
fuselage, wing, control  surfaces  and  tip  missile  configuration  is  analyzed.  Neither  Ref.  10  (used  for 
constructing  the  analytical  wing  model)  nor  Ref.  54  specify  whether  the  wing  is  empty  or  full  with 
fuel.  A  difference  in  wing  fuel  level  might  explain  the  slight  difference  in  first  wing  bending 
frequency  between  analysis  and  test. 
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EQUIVALENT  PLATE  CAPABILITY  TESTING 


NATURAL  FREQUENCIES  ( HZ )  OF  THE  YF16 


CANTILEVERED  WING/LE  FLAP / 
FLAPERON  ASSEMBLY 


F-F  A/C  WITH  WING  TIP 
MISSILES  (ANTI-SYMMETRIC) 


No. 

F.E.M 

TSO 

PRESENT 

No. 

GVT 

PRESENT 

(REF.  10) 

(REF.  10) 

CODE 

(REF. 54) 

CODE 

1 

10.67 

10.74 

9.98 

1 

6.5 

6.39 

2 

33.92 

35.05 

34.98 

(Missile 

3 

35.78 

42.75 

36.48 

Pitch) 

4 

56.45 

64.24 

54.02 

2 

8.0 

8.70 

5 

62.47 

73.43 

65.28 

(Wing  1st 

6 

67.96 

95.31 

73.57 

Bending) 

Figure  IQ 


-  11  - 


NUMERICAL  TESTING  (CONTINUED1) 


The  first  six  mode  shapes  for  the  cantilevered  YF16  example  (without  tip  missile),  generated  by 
the  new  multi-element  equivalent  plate  analysis,  are  shown  in  Fig.  11.  These  mode  shapes 
correlate  well  with  finite  element  results  reported  in  Ref.  10.  The  quality  of  this  correlation  can  be 
attributed  to  the  high  order  of  control  surface  displacement  representation  and  better  modeling  of 
elastic  point  attachment  of  the  control  surfaces  to  the  wing  box. 

A  typical  computation  time  for  the  static  analysis  of  the  wing  of  Ref.  43  (including  the  calculation 
of  384  displacement,  slope  and  stress  constraints  and  their  sensitivities  with  respect  to  inner  and 
outer  panel  skin  thicknesses  at  an  array  of  points  over  the  wing)  is  12.6  epu  seconds  on  the 
UCLA  IBM  3090.  Analysis  and  constraint  generation  for  YF16  six  static  load  cases  and  natural 
modes  takes  18.9  seconds.  These  relatively  short  computation  times  are  essential  to  the 
construction  of  an  efficient  multidisciplinary  synthesis  capability. 
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MODE  SHAPES  OF  THE  YF16  CANTILEVERED  WING 


MODE  5  MODE  6 


Figure  11 

YL 


LIFTING  SURFACE  UNSTEADY  AERODYNAMICS  :THE  PCKF  METHOD 


Along  the  line  of  improving  the  mathematical  modeling  of  the  scrvoaeroclastic  wing  dynamic 
system,  the  use  of  lifting  surface  theory  (Ref.  46,47)  for  the  calculation  of  the  unsteady 
aerodynamic  loads  is  considered  a  definite  step  forward  compared  with  strip  theories.  Lifting 
surface  aerodynamics  are  still  widely  accepted  in  the  aerospace  industry  for  the  flutter  and  gust 
response  analysis  of  airplanes  in  the  subsonic  and  supersonic  speed  regimes.  Thus  including  lifting 
surface  modeling  in  the  analysis  part  of  a  multidisciplinary  wing  synthesis  is  important  if  the 
synthesis  of  real  wings  is  sought. 

In  the  PCKF  method  for  the  solution  of  the  integral  equation  relating  downwash  and  pressure 
distribution  over  a  lifting  surface  (Ref.  48-51)  an  assembly  of  lifting  surfaces  is  divided  into  a 
group  of  trapezoidal  boxes,  as  shown  in  Fig.  12  for  a  subsonic  case. 
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y£nNG_SURjEACE_UNSTEADY 


.nonnYNAMlCSjTHgJC^lglHQD 


iTInN  BY  AN  ASSEMBLY  OF  TRAPEZOIDAL  BOXES 
MODELING  a  CONFIGLRATION  B 
(SUBSONIC) 


THE  PCKE  METHOD  :  SOME  ANALYTICAL  ASPECTS 


The  pressure  distribution  on  each  box  is  approximated  by  weighting  functions  representing  the 
known  pressure  singularities  along  the  box  edges  multiplied  by  a  series  of  polynomials  orthogonal 
to  these  weighting  functions.  Collocation  points  over  the  planform  are  chosen  so  as  to  minimize 
the  error  in  the  pressure  integrals  needed  to  calculate  generalized  aerodynamic  forces.  The  PCKE 
method  is  fast,  accurate  and  especially  suited  to  handle  wing/control  surface  configurations.  It  is 
more  accurate  than  the  vortex  lattice  method  especially  when  leading  edge  flaps  or  controls  with 
gaps  around  them  are  considered  (Ref.  48).  This  is  due  to  the  inability  of  lattice  methods  to 
impose  the  pressure  singularities  along  the  different  boundaries  of  the  wing.  In  the  present 
application  it  is  integrated  with  the  equivalent  plate  structural  analysis  to  generate  a  set  of 
generalized  loads  for  the  same  generalized  polynomial  coordinates  used  for  structural  analysis.  The 
number  of  collocation  points  per  box  and  the  number  of  integration  points  used  are  carefully 
selected  to  be  compatible  with  the  order  of  displacement  polynomials  used. 


THE  PCKF  METHOD  :  ANALYTICAL  ASPECTS 


POLYNOMIAL  SERIES  APPROXIMATION  FOR  PRESSURE  OVER  A  BOX  : 


Mv>l) 


ns  pan 


nctotd 


i- I 


impJLn)  MZ)p,U) 


where  m  =  (j—  \)nchord  +  i 


WEIGHT  FUNCTIONS  REPRESENTING  KNOWN  PRESSURE  SINGULARITY  ALONG 
BOX  EDGES 


COLLOCATION  POINT  PLACEMENT  :  OPTIMAL  SO  AS  TO  MINIMIZE  ERROR  IN 
PRESSURE  INTEGRALS  (  GENERALIZED  AERODYNAMIC  FORCES  ) 


ADVANTAGES  : 


SUBSONIC  /  SUPERSONIC 

GENERAL  NON  PLANAR  WING/CONTROL  SURFACE  CONFIGURATIONS 

FAST  CONVERGENCE  OF  GENERALIZED  LOADS  WITH  INCREASED  NUMBER  OF 
POLYNOMIALS  -  HIGH  COMPUTATIONAL  SPEED 

GOOD  ACCURACY  OF  CONTROL  SURFACE  HINGE  MOMENT  AND  CONTROL 
SURFACE  DERIVATIVES  (VORTEX  LATTICE  METHOD  OVERPREDICTS  HINGE 
MOMENTS)  -  IMPORTANT  FOR  SERVOAEROELASTIC  MODELING 

EXTENSIVE  NUMERICAL  TESTING  BY  THE  DEVELOPERS  FOLLOWED  BY 
ACCURATE  RESULTS  IN  THE  FLUTTER  AND  SERVOAEROELASTIC  ANALYSIS  OF 
THE  F16  IN  AN  INDUSTRY  ENVIRONMENT 

A  DEFINITE  IMPROVEMENT  OVER  STRIP  THEORIES 


Figure  13 
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UNSTEADY  AERODYNAMICS  FINITE  STATE  MODELING 


The  generalized  aerodynamic  loads  in  the  Laplace  transformed  small  perturbation  equations  of 
motion  (for  a  steady  level  flight)  given  below  are  transcendental  functions  of  the  Laplace  variable 
s.  The  flight  dynamic  pressure  and  flight  speed  are  qD,  respectively;  M,C,K  are  the 
mass, damping  and  stiffness  matrices:  Q(s.\1ach)  is  the  matrix  of  generalized  aerodynamic  forces  in 
the  Laplace  domain:  is  the  Laplace  transformed  vertical  gust  velocity  ;  S  is  a  reference  area 

and  is  the  vector  of  Laplace  transformed  generalized  displacements. 

To  use  modem  control  system  analysis  and  design  techniques,  it  is  necessary  to  cast  them  in 
Linear  Time  Invariant  (L'l'I)  state  space  form.  The  common  practice  is  to  match  rational  function 
approximations  to  generalized  aerodynamic  loads  calculated  for  harmonic  motion  at  a  set  of 
reduced  frequencies  (Ref.  53).  There  is  a  resulting  increase  in  the  order  of  the  LTI  state  space 
model  due  to  the  addition  of  aerodynamic  states.  This  increase  in  size  can  be  quite  significant. 
With  n  generalized  displacements,  each  lag  term  in  the  commonly  used  Roger  approximation  (see 
Ref.  53  for  further  detail)  adds  n  states  to  the  model  order.  Since  four  lag  terms  are  usually  needed 
for  a  reasonable  approximation  in  this  method,  4n  states  are  added  to  the  system.  This  makes  it 
computationally  expensive  to  carry  out  any  control  system  analysis  and  behavior  sensitivity 
analysis  using  state  space  techniques.  In  the  Minimum  State  Method  of  Karpel  (Ref.  52  ),  the 
functional  dependence  of  the  generalized  aerodynamic  force  matrix  on  the  Laplace  variable,  is 
approximated  by  a  rational  expression  of  a  special  form  so  as  to  reduce  the  number  of  added 
states  needed  to  achieve  given  quality  of  fit. 

Given  the  generalized  aero  forces  in  simple  harmonic  motion  for  a  number  of  reduced  frequencies, 
it  is  possible  to  match  the  approximation  exactly  to  the  data  for  k  =  0  and  one  other  reduced 
frequency.  This  determines  the  matrices  P„  P2,  P3.  Chossing  R  to  be  a  diagonal  matrix  with 
negative  elements,  the  matrices  D  and  E  arc  determined  in  an  iterative  process  so  that  the 
approximation  fits  the  rest  of  the  data  in  a  least  square  manner. 
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UNSTEADY  AERODYNAMICS  FINITE  STATE  MODELING 


THE  SMALL  PERTURBATIONS  LAPLACE  TRANSFORMED  EQ.  OF  MOTION  OF  AN 
ELASTIC  AIRPLANE  IN  LEVEL  FLIGHT  : 

{[A/]*2  +  [C]j+  [/Q}{<j(5)}  -  qD S  [QLsMamm)  =  Id*  {Qds,Mach))^- 


PURPOSE  OF  FINITE  STATE  MODELING  :  CAST  EQ.  OF  MOTION  IN  LINEAR 
TIME  INVARIANT  STATE  SPACE  FORM 

PRINCIPLE  :  RATIONAL  FUNCTION  APPROXIMATIONS  OF  UNSTEADY 

AIRLOADS  IN  TERMS  OF  LAPLACE  VARIABLE 

PRICE  :  ADDED  STATES 

MINIMUM  STATE  APPROXIMATION  FORM  : 


eg,)]  =  C py  +  iP{\s  +  [P3]  +  [D][ii  -  *]-'[£]* 

MATCHING  PROCESS 

♦  GENERALIZED  AERO  FORCES  ARE  GIVEN  FOR  HARMONIC  MOTION  AT  A  SET 
OF  REDUCED  FREQUENCIES 

♦  A  SET  OF  AERODYNAMIC  LAG  TERMS  IS  CHOSEN  :  R,t 

♦  P,  IS  EQUATED  TO  0  k  =  0  ) 

♦  P„  P2  ARE  EXPRESSED  IN  TERMS  OF  D,E  SO  AS  TO  ENSURE  PERFECT  FIT  AT  A 
SELECTED  REDUCED  FREQUENCY  k, 

♦  D,E  ARE  DETERMINED  IN  AN  ITERATIVE  LEAST  SQUARE  PROCESS  TO  FIT 
THE  REST  OF  THE  DATA 


ADVANTAGE : 

MINIMAL  INCREASE  IN  MODEL  ORDER 
PROBLEMS  : 

ITERATIVE  PROCESS  IS  TIME  CONSUMING 

RELATIVELY  LITTLE  EXPERIENCE  WITH  REAL  CONFIGURATIONS 


Figure  14 


SOME  PRELIMINARY  MINIMUM  STATE  FITS  TO  A  LARGE  MATRIX 
OF  UNSTEADY  AERODYNAMIC  GENERALIZED  FORCES 


Order  reduction  of  the  state  space  model  used  for  servoaeroelastic  stability  and  control  analysis  is 
essential  for  synthesis  purposes  in  order  to  make  the  analysis  cycle  as  computationally  fast  as 
possible.  This  motivates  the  choice  of  the  Minimum  State  Approximation  for  finite  state  unsteady 
aerodynamic  modeling  in  the  current  research.  Preliminary  tests  of  the  quality  of  approximation 
achieved  when  applied  to  a  large  matrix  of  generalized  aerodynamic  forces  show  promising  results. 
A  44  x  44  matrix  of  generalized  aero  forces  for  the  YF16  airplane  with  tip  missiles  is 
approximated  using  only  22  lag  terms.  Comparison  with  a  one  lag  term  Roger  approximation 
(which  will  add  44  aerodynamic  states  to  the  model)  shows  an  advantage  of  the  minimum  state 
approach. 
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SOME  PRELIMINARY  MIN  STATE  FITS  TO  A  LARGE  MATRIX 
"(TF  UNSTEADY  aER6pynami£  GENERALIZED  FORCES' 


SOME  RECENT  EXAMPLES  OF  QUALITY  OF  FIT  FOR  A  YFI6  COMPLETE  A/C 
CONFIGURATION  : 

(  44  POLYNOMIAL  GENERALIZED  COORDINATES  ) 

SOME  LOW  ORDER  SHAPE  FUNCTIONS  FOR  THE  WING  BOX  : 
f  (*j0«  L  A  (jf«v) -  x,  A  (*«?)« ** 

A  HIGHER  ORDER  SHAPE  FUNCTION  :fu  (xjO-  * 


ROGER  APPROXIMATION  BASED  ON  1  LAG  (44  ADDED  STATES) 
MINIMUM  STATE  BASED  ON  22  LAGS  (22  ADDED  STATES) 


I  c 


Alii: 


CONTROL  SYSTEM  MODELING 


A  block  diagram  of  the  actively  controlled  scrvoaeroelastic  system  is  shown  in  figure  16.  Airplane 
motions  (acceleration  and  angular  rates)  are  sensed  by  a  set  of  sensors  placed  at  different  points 
on  the  structure.  The  resulting  signals  are  used  as  inputs  to  the  control  law  block  which 
commands  control  surface  actuators.  The  control  surface  motions  guarantee  stability  and 
desirable  dynamic  response  of  the  complete  system. 

For  the  control  system,  only  sizing  type  design  variables  are  considered  at  present  to  keep  the 
balance  in  our  approach,  and  these  are  the  coefficients  of  numerator  and  denominator 
polynomials  in  the  control  law  transfer  functions.  Control  surface  locations,  sensor  locations,  the 
structure  of  the  control  system  and  order  of  transfer  functions  are  preassigned.  It  is  assumed  that 
sensor  and  actuator  transfer  function  are  given,  although  the  formulation  is  general  enough  to 
allow  treating  their  elements  as  design  variables  as  well. 
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CONTROL  SYSTEM  MODELING 


gust  excitation 


LTI  STATE  SPACE  MODEL  AND  STABILITY  ANALYSIS 


Formulations  of  the  state  space  control  augmented  servoacroeiastic  equations  of  motions  can  be 
found  in  many  works  on  active  flutter  suppression  (e.g.  Ref.  55).  A  transfer  function  model  of  an 
element  of  the  control  system  (whether  sensor,  actuator  or  a  control  law)  can  be  transformed  into 
a  state  space  model,  where  the  A,  fi,  C  and  D  matrices  are  explicitly  expressed  in  terms  of  the 
transfer  function  numerator  and  denominator  polynomial  coefficients.  Assembly  of  the  sensor, 
actuator,  control  law,  structural  dynamics,  gust  and  unsteady  aerodynamics  state  space  models 
leads  to  the  system  matrices  U.  V  and  W  in  a  LTI  state  space  model  of  the  whole  system.  These 
matrices  are  functions  of  the  structural  design  variables  through  their  dependence  on  the  stiffness 
and  mass  matrices.  They  depend  on  the  control  system  design  variables  through  their  dependence 
on  the  state  space  models  of  the  control  elements. 

For  given  flight  conditions  (Mach  number  and  Altitude)  the  stability  of  the  system  is  determined 
by  the  real  part  of  the  eigenvalues  of  a  generalized  eigenvalue  problem.  Sensitivity  of  a  critical 
eigenvalue  with  respect  to  any  design  variable,  p  ,  is  calculated  using  standard  eigenvalue 
sensitivity  analysis  based  on  the  derivatives  :  cUjdp,  dVfdp  and  the  left  and  right  corresponding 
eigenvectors  {rj/},  {<p}.  It  is  planned  to  use  the  original  Ritz  functions  directly  as  generalized 
coordinates.  This  approach  leads  to  an  increased  order  model  but  avoids  natural  mode 
calculation  and  aerodynamic  force  updates  associated  with  natural  mode  reduced  models. 
Computation  times  and  accuracy  will  determine  whether  there  is  a  need  to  resort  to  natural 
modes.  Alternative  approximations  to  system  eigenvalues  in  terms  of  structural  and  control 
system  design  variables  will  be  studied. 
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LTI  STATE  SPACE  MODELING  AND  STABILITY  ANALYSIS 


STATE  SPACE  MODELS  OF  ACTUATORS,  SENSORS  AND  GUST  FILTER  : 


s  l*/}  =  C  4  3  to }  +  lBH  {q} 


O'/}  =  C  c,  ]  {x,} 

i=  ACT  FOR  ACTUATORS 
i=S£,V  FOR  SENSORS 
/'  =  C  FOR  GUST 

STATE  SPACE  MODEL  OF  THE  CONTROL  BLOCK  : 


s  (xLAl v)  ~  t  3  ixLAw)  +  C BLA\v]  iuLAw) 


iyiAw)  ~  [  Claw]{xLAw)  +  C ^LAwl{uLAw) 

THE  A,B,C,D  MATRICES  ARE  EXPLICITLY  EXPRESSED  AS  A  FUNCTION  OF  THE 
CONTROL  SYSTEM  DESIGN  VARIABLES. 


SYSTEM  STATE  VECTOR  {  x  }  =  {x,  x2  ..  x,  }  CONTAINS  : 

STRUCTURAL  STATES  ;  ACTUATOR  STATES  ;  SENSOR  STATES  ;  CONTROL  LAW 
STATES  ;  GUST  STATES  ;  AERODYNAMIC  STATES  ASSOCIATED  WITH 
GENERALIZED  AERO  MATRIX  ;  AERODYNAMIC  STATES  ASSOCIATED  WITH 
GUST  VECTOR  . 


THE  CLOSED  LOOP  STATE  SPACE  EQUATIONS  OF  THE  COMPLETE  SYSTEM  : 
sm  =  tnwi)}  +  mu^s) 


STABILITY  BY  EIGENVALUE  ANALYSIS  : 

A  [£/(p)]  {tf}=  LV(p)U4>) 


EIGENVALUE  SENSITIVITY  WITH  RESPECT  TO  DESIGN  VARIABLE  p  : 


ii 

dp 


« l/T  U<(> 


Figure  17 
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STATUS  OF  MULTIDISCIPLINARY  ANALYSIS  AND 
BEHAVIOR  SENSITIVITY 


Figure  1 8  presents  status  of  research  activities  associated  with  the  development  of  the  analysis  and 
sensitivity  capabilities  for  the  multidisciplinary'  synthesis  of  wings.  It  is  expected  that  based  on 
these  capabilities,  it  will  be  practical  to  synthesize  on  a  preliminary  design  level  realistic 
representations  of  control  augmented  wings.  The  generality  of  the  approximation  concepts  based 
mathematical  programming  aproach  to  synthesis  and  the  realism  in  modeling  are  expected  to  be 
of  major  importance  in  coping  with  complicated  multidisciplinary  interaction,  where  little 
experience  exists  and  intuition  is  often  misleading. 


ANALYSIS  AND  BEHAVIOR  SENSITIVITY 
STATUS 


STRUCTURE  AERODYNAMICS  *  CONTROL 


FORMULATION  :  ANALYSIS  + 

SENSITIVITY  + 

ANALYSIS  IMPLEMENTATION  :  + 

ANALYSIS  TESTING  :  + 

SENSITIVITY  IMPLEMENTATION  :  + 

SENSITIVITY  TESTING  :  + 

APPROXIMATION  CONCEPTS 
PERFORMANCE  ASSESSED  :  + 


+ 

+ 

+ 

+ 

in  progress 


+ 

+ 

+ 

in  progress 
in  progress 


♦  AERODYNAMICS  INCLUDE  : 

♦  UNSTEADY  AERODYNAMICS  FOR  SERVOAEROELASTIC  ANALYSIS 

♦  STEADY  TRIM  AND  DRAG  CALCULATIONS 


Figure  18 
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